################################################################################
# PACKAGE AND DATA LOADING
################################################################################

pacman::p_load(sf,rnaturalearth,lwgeom,tmap)

borders <- st_read('land_borders.shp')
walls <- st_read('europe_walls.shp')

world <- ne_countries(scale = "medium", returnclass = "sf")
world = lwgeom::st_transform_proj(world, crs = "+proj=moll")
bbox_new <- st_bbox(c(xmax = 3530000, xmin = (-1080000), ymax = 7960000, ymin = 4130000), 
                    crs = st_crs("+proj=moll"))
xrange <- bbox_new$xmax - bbox_new$xmin
yrange <- bbox_new$ymax - bbox_new$ymin
bbox_new <- bbox_new %>%  
  st_as_sfc()

################################################################################
# FIGURE 4
################################################################################

pdf("Figure 4.pdf",750,1125,res=150)
tm_shape(world,bbox=bbox_new) +
  tm_polygons(col="black",border.col="black") +
  tm_shape(borders) + tm_lines("steelblue4",lwd=.5) +
  tm_shape(walls) + tm_lines("darkorange",lwd=3.5)
dev.off()
